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A theory is presented for Ac, the coefficient of the first-order correction in the density of the 
collective diffusion coefficient, for protein spheres interacting by electrostatic and adhesive forces. 
An extensive numerical analysis of the Stokesian hydrodynamics of two moving spheres is given 
so as to gauge the precise impact of lubrication forces. An effective stickiness is introduced and a 
simple formula for Ac in terms of this variable is put forward. A precise though more elaborate ap- 
proximation for Ac is also developed. These and numerically exact expressions for Ac are compared 
with experimental data on lysozyme at pH 4.5 and a range of ionic strengths between 0.05 M and 
2 M. 



I. INTRODUCTION 

Pick's first law states that the particle flux is equal to 
minus the collective diffusion coefficient times the gradi- 
ent of the particle concentration. For colloids or macro- 
molecules in solution, this collective (also called cooper- 
ative or mutual) diffusion coefficient is often determined 
experimentally with the help of dynamic light scattering. 
If one extrapolates this coefficient to a vanishing concen- 
tration of particles, it reduces to the single-particle diffu- 
sion coefficient since the interactions between the parti- 
cles are presumably negligible then. At non-zero volume 
fractions, particle interactions, such as those of electro- 
static and hydrodynamic origin, influence the diffusion. 
At low enough concentrations, where three- and higher 
body interactions may be disregarded, the parameter Ac 
characterizes the departure from the single-particle re- 
sult. 

The concentration dependence of the collec- 
tive diffusion coefficient of proteins has been 
studied extensively in experiments, for example 
in the case of hemoglobirtii^i^iiii^, bovine serum 
albumin^iii^i^, /3-lactoglobuliniS, ovalbumirtii and 
lysozymei^^'i^^'il^'ia^^^. On the theoretical 
side, a fair number of papers2^i^i2^i^i21i^i2Si^i2i deal 
with the diffusion of interacting colloidal particles in 
solution. Apart from giving insight into the diffusion 
as such, the coefficient Ac is also important because 
it could yield information about the complex pair 
interaction between protein molecules. Moreover, it has 
been argued that Ac may be an alternative parameter 
useful in diagnosing under what conditions proteins 
would crystallizeii. 

In Ref. [s^ we approximated globular proteins in wa- 
ter with added monovalent salt by hard spherical parti- 
cles that interact through a short-range attraction and 
a screened electrostatic repulsion. We appropriately re- 



1 Address for correspondence: T. Odijk, P.O. Box 11036, 2301 EA 
Leiden, the Netherlands. E-mail: odijktcf@wanadoo.nl 



placed this system by one of spherical particles with 
sticky interactions only. At infinite dilution the effective 
stickiness is readily determined by equating the respec- 
tive second virial coefficients of the two systems. In the 
effective stickiness, part of the bare adhesion is balanced 
against the electrostatic repulsion. 

In the next section, we formulate a theory for the coef- 
ficient Ac- We first introduce the interaction used previ- 
ously to compute protein solution properties^ and give 
expressions for the effective stickiness. We then outline 
the formal expression for Ac due to Pelderho:^ in terms 
of the pair potential between two protein spheres and 
a hydrodynamic mobility function. Although the latter 
has been studied often in the past, we present a more ex- 
tensive numerical analysis in order to gain more insight 
into the asymptotics of the lubrication regime for two 
moving spheres very close to each other. The coefficient 
Ac is then computed in three ways: exactly via numer- 
ics and in terms of two convenient approximations. In 
section III, we compare these predictions for Ac with ex- 
periment. A discussion of the results is given in the last 
section. 



II. THEORY 

A. Effective interaction 

We model the globular proteins as spherical particles 
of radius a with a total charge Zq per particle that is 
uniformly distributed over its surface. Here q is the el- 
ementary (proton) charge. Por convenience, we scale all 
distances by the radius a and all energies by ksT where 
/cs is Boltzmann's constant and T is the temperature. 
We approximate the interaction between two proteins by 
a steric repulsion plus a short-range attraction of scaled 
range and constant absolute magnitude Uai and a 

far-field Dcbye-Hiickel potential. The latter describes the 
Coulomb repulsion that is screened due to the presence of 
monovalent salt of ionic strength /. The effective number 
Zef f of charges associated with the far field is computed 
in the Poisson-Boltzmann approximation. See Refs. [s^ 
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and I33I for further details. The total interaction Ut{x) 
between the two particles with center-of-mass separation 
r is thus of the form 

r 00 < a; < 2 

Urix) ^ { Udh{x)-Ua 2<x<2 + 5 , (1) 

i UoHix) x>2 + 6 



Here, the Dcbyc-Hiickel interaction is given by 



UDHix)^2^- 



(2) 



(3) 



where ^ = ^ (^"y+^) fj. = na. The Debye length 

is defined by = SttQI and the Bjerrum length 
by Q = (f /eksT, which equals 0.71 nm in water at 298 
K (e is the permittivity of water); ^ = 3.28a-\/j, if the 
radius a is given in nm and the ionic strength I in M. We 
suppose 1-1 electrolyte has been added in excess so / is 
the concentration of added salt only. We have derived a 
pcrturbativc expression for the effective charge qZf.f f in 
the Poisson-Boltzmann approximation^ 



/i 



c3'^£;i(3m). (4) 



Here, Ei{x) is the exponential integral defined by 
Ei{x) = dtt~^e~* . Eq. (HI is numerically consis- 
tent with a different form recently proposed by Aubouy 
et al<^ which is also valid at large values of Z . 

We want to replace the system of particles interact- 
ing through the complicated interaction ([T]) by a system 
of particles interacting through a simpler potential, the 
adhesive hard sphere (AHS) potential of Baxter— 
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< a; < 2 
2 < a; < 2 
x>2 + uj 



(5) 



Here, r is a positive constant which signifies the strength 
of the effective adhesion and the limit a; J, has to be 
taken appropriately after formal integrations. In order 
to replace the original system by this simpler system, we 
have to find the correspondence between the parameter 
T in the AHS potential and the parameters ^, fi, 6 and 
Ua in the original interaction Eq. ([T]). In this case, we 
do this by matching the respective second virial coeffi- 
cients, which ensures that the free energy of the two sys- 
tems at small concentrations are identical. We emphasize 
that in the general case, at arbitrary concentrations, we 
have to match the complete free energies of the respective 
systems^i^^; it is then incorrect to focus on the second 
virials as has often been done in the past. 



B. Stickiness parameter 



We already determined the stickiness parameter r in 
a previous paper—. Here we reproduce the main results. 



The second virial coefficient B2 is given by 
B, = - / dr f 1 - e-^( 



(6) 



where U (r) is the pair potential scaled by kgT, and r 
is the unsealed position vector connecting the centers of 
mass of the two particles. For the pair interaction of 
Eq. ID), B2 may be expressed by 



B2 = Bf ^ ( 1 + ^ J 



(7) 



where we introduce the following integrals 
J = dxx^ (1 - e-^'-^^)) = Ji ~ (e^^* - 1) J2, (8) 



Ji 



dxx^ (l_e"'^°«("=) 
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Here, — 167ra'^/3 is the value of B2 if the pro- 

teins were solely hard spheres and a ~ - — We 
equate Eq. ^ with the second virial coefficient of the 
AHS model 



B2 = B2 



HS 



1 



which results in a stickiness parameter r given by 

2 



(11) 



(12) 



From Eqs. ([T]) and ([8]) we see how part of the original 
attraction is compensated by repulsive electrostatics. 



C. General expression for Ac 

For small volume fractions (p of spherical particles, the 
collective diffusion coefficient Dc may be written as 



Dc = Do{\ + \c^ + 0{4>'')) 



(13) 



where Dq is the diffusion coefficient in the dilute limit. 
The linear coefficient Ac may be split up into five 
contributions^^ 



Ac — Ay + Ao + Ad + A5 + A^. 



(14) 
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These terms have been studied for some time^^i^iiS^iS^: 
there is a virial correction because a fluctuation in the 
osmotic pressure drives difi^usion 



Ay = 3 / dxx" [1 — c 



-U{x] 



(15) 



and four terms arising from the mutual friction between 
two hydrodynamically interacting spheres. An Oseen 
contribution 



Ao = 3 / da;a;(e-^(^) - l) , 



and a dipolar contribution 

\d = 1, 



(16) 



(17) 



express the long-range hydrodynamic interaction be- 
tween two particles 1 and 2 whereas the short-range part 
of the hydrodynamic interaction comes into play in the 
term 



3 

X] 
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Finally, the modification of the single-particle mobility is 
expressed by 

da; x2c-^(^) {A[\ (x) + 2Bl\ (x)) . (19) 

Here, A\\{x), A\*2ix), B{\{x) and Bfi^ix) arc dimension- 
less hydrodynamic functions given in terms of the trans- 
lational mobility matrix for two spheres centered at Ri 
and i?2 {r ^ Ri — R2) and acquiring velocities Vi and V2 
as a result of the forces Fi and F2 acting on the spheres 

V, ^ fil\{l,2) ■ F, + fi{\{l,2) ■ F2 (20) 

V2 ^ t^l\il,2) ■ Fi + ^ll\il,2) ■ F2 (21) 
In the notation of Cichocki and Felderho^, we have 
1 



Mn(l,2) 



M*i2(l,2) 



(22) 
(23) 



where 77 is the viscosity of the solvent and / is the unit 
tensor. The mobility tensors in Eq. (PT|) are given by in- 
terchanging the labels in Eqs. and while taking 
into account the symmetry relations 



(24) 



Recall that the particles have a hard-core interaction 
for a; < 2 so exp —U (x) vanishes for x < 2. We then sum 
Eqs. p^ - fTO]) and conveniently rewrite Ac as follows 



Ac = Co + ci 



dx : 



(l-e-^(-)) + 



R. 



(25) 



The constant cq equals the value Ac would adopt if the 
spheres were hard but without any other interaction 

Co = 3 / dx x^ — 3 f dxx + 1 + f dxx^ ( h{x) 

^0 Jo J2 V ^ 

3 + ^ dxx^ (^h{x) - . (26) 



Here, h{x) is the sum of scalar mobility functions 

h{x) = A\\{x) + A\*2{x) + 2Bfi (x) + 25^2(2;). (27) 



The residual term R in Eq. ((251) depends on the actual 
interaction 

R= /"°°dxa;2 (e~^(") - 1) (/^(.t) -/i(2)) (28) 



though it would vanish if the interaction U were adhesive 
and purely of the Baxter type (see Eq. ([5])). The second 
term on the right hand side of Eq. (|25p is proportional 
to the constant 



ci = 3 - h{2) 



(29) 



and the integral is related to the second virial coefficient 
B2 by (see Eqs. d?]) and ©) 



dxx^ l-e-^'^-A ^- 



f B2 



' 3\B^' 
The resulting expression for Ac is 
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Ac = Co + — 



8ci / B; 



3 \B^s 



-1]+R 



(30) 



(31) 



which we can evaluate once we know h{x) given by 
Eq. (EH). 



D. Hydrodynamics 

The function h{x) was discussed by Batchelor— in 
his theory of the diffusion of hard spheres. The sum 
Afi + A^i2 pertains to the mobility of a pair of spheres 
moving in the direction of their line of centers whereas 
Bu 

perpendicular to that line 
itt _L 1 n. . = Rtt 



B*2 is related to their mobility when they move 



(Note that in Rcf. |3 
itt 



All = 

Al\ + 1, Bn = Bl{ + 1, A12 = AY^ and B^ = B[l). In 
the latter case, because the spheres are couple-free, the 
spheres must rotate as the pair translates. At small sepa- 
rations {x — 2 <C 1), lubrication forces with a logarithmic 
singularity ln~^(a; — 2) are then expected to develop on 
general grounds^. Goldman et al^ proposed a form for 
the singularity which we will test below. 

Batchelor— computed h{2) = 1.312 on the basis of nu- 
merical work on the mobilities of touching spheres^SilL. 
Cichocki and FeldcrhoS^ evaluated co = 1.454 (Eq. (gB])) 
by numerically summing their series expansions of the 
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hydrodynamic interactions while keeping track of a log- 
arithmic singularity at close separations. Here we reana- 
lyze h{x) and go well beyond previous computations^^i^ 
in order to gain more insight into the nature of the sin- 
gularity and to calculate the residual R. 

We assume the interaction U (x) is of short range so we 
focus only on h{x) for x — 2 < 1. First, we get an expres- 



the following analytical expression for h{x) valid at small 
separations 



sion for A\\ 



A^i2 an infinite sum from the results of 



Stimson and Jeffery^ who expressed the hydrodynamic 
problem in terms of bispherical coordinates. (Note that 
there is an error in their paper as pointed out in, for ex- 
ample, Ref-IUin which one may find a similar expression 
for A^ii — A*i2 in case one needs and separately) . 
Calculating Bl\ and is more involved. We use the 
numerical scheme by O'Neill and Majumdar— which is 
similar to that of Goldman ct al^. (Note that there are 
a few typographical errors in Ref. [4^. In their Eq. (3.9) 
d should be di, the expression for v in Eq. (4.1) should 
have a minus sign, ^, and ip in Eqs. (4.3)-(4.5) should 
be replaced by c^, ccj) and c-0 respectively, and sinh^|/3| in 
Eq. (5.10) shouldbesinh^l^l. Also, to obtain i:»2(A„, S„) 
(Eq. (3.29)) from Di{An,Bn) (Eq. (3.28)) the signs of 
6n-i, Sn and Sn+i should be reversed as well (we only 
checked the case of spheres of equal size). Their Table I 
is correct, however, for spheres of equal size, apart from 
the value for (712(1, 0.1) which should read -0.1017 instead 
of -1.1017). 

In order to investigate the regime of lubrication for a 
pair of spheres moving under the action of applied forces 
normal to their line of centers, we performed the numer- 
ical analysis down to r/a — 2 ~ 10~^° which implies two 
million terms in the series expansions are needed. We at- 
tempted to speed up the iteration by adapting the recur- 
rence relationships introduced more recently by O'Neill 
and Bhatt^ for a sphere moving near a wall to the case 
of two spheres. However, this did not turn out to be 
useful as it is for the wall configuration^. One way of 
circumventing series expansions could be to elaborate on 
the trial functions initially used by Fixman in his varia- 
tional theorem for the mobility matrix-^ but we did not 
investigate this. 

Goldman et al^ were the first to give a comprehensive 
analysis of the mobility of a pair of identical spheres of ar- 
bitrary orientation. They numerically solved the Stokes 
and continuity equations using expansions in terms of 
bipolar coordinates to high order. For moving spheres 
whose line of centers is perpendicular to the applied force, 
the force consists not only of a term arising from pure 
translation but also a term stemming from pure rotation 
of the spheres. The latter involves a torque on one sphere 
diverging aa^ 



Tr 



3 ln(.T - 2) 



(32) 



at very small separations where Q is its angular veloc- 
ity. Eq. (|32p was derived by extending the nontrivial 
lubrication theory of Ref. |49| in which inner and outer 
regions have to be matched. Eq. (|32p ultimately leads to 



h{x) = h{2) 



0.47666 



\u{x - 2) + C2 



+ 0{x-2). 



(33) 



The coefficient 0.47666 is computed from the numerical 
tables presented in Ref. |3^. We have added a constant 
C2 to the logarithm because we expect the next higher 
order term in Eq. (j32[) to be a constant judging by the 
earlier analysis of the sphere-wall problem^. In Fig. [T] 
we have fitted Eq. (|33|1 to the numerical results discussed 
above, letting h{2) and C2 be adjustable. The intercept 
h{2) = 1.30993 turns out to be close to the value 1.312 
quoted above for touching spheres which lends credence 
to the validity of the asymptotic expression that we pro- 
pose. Moreover, the resulting coefficient C2 = —4.694 and 
the concomitant shift in Eq. ([32)1 are consistent with the 
numerical values of the torque at small separations as 
presented in table 3 of Ref. [S^. 
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FIG. 1: The hydrodynamic function h plotted in terms of 
the variable s = — l/(ln(x- — 2) — 4.694). Squares denote 
results from the numerical analysis to the accuracy as ex- 
plained in the text. The straight line signifies the function 
h = 0.47666s + 1.30993. 



Next, we derive an expression for the residual term 
given by Eq. (|28p . First, we propose an initial estimate 
/io(a;) for h{x). We have plotted the numerical values 
of h{x) as a function of x in Fig. [21 As a result of the 
lubrication regime, h has a maximum as displayed in the 
inset. However, h{x) is only a strongly varying function 
for X < 2.04. We therefore simply force a linear fit to the 
data for h at x = 2.1, 2.2 and 2.3 



ha{x) « 1.3670 - 0.4745(x - 2). 



(34) 

We then insert this estimate into Eq. (pS)) and add a cor- 
rection term so as to derive an expression for R accurate 
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enough for our purposes. 



-0.147 



/ B2 



1 



0.4745 J dxx^{x-2)(^ 
9x 10-^ (l_e-^(2)j^ 



1 



-U{x] 



(35) 



The first term on the right comes from the fact that the 
linear interpolation gives ho{2) = 1.3670 whereas the real 
value is h{2) = 1.312. Since the interaction usually does 
not change appreciably for 2 < a; < 2.04, it is straightfor- 
ward to write an estimate for the error — the third term — 
owing to the deviation of Eq. (pij) from the exact function 
h{x) (see inset Fig. In our case the error term turns 
out to be an order of magnitude smaller than the first 
two terms. 
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FIG. 2: The hydrodynamic function h{x) as a function of the 
dimensionless separation x = r/a between the centers of two 
spheres. The straight line signifies ho{x) given by Eq. (|34|l . 



E. Determination of Ac 



It is clear from Eq. (|28p that R would vanish if the 
actual interaction were a pure AHS potential. If we then 
insert Eq. ^ into Eq. we obtaiii^^ 



(36) 



, 2ci 
Ac = c„-— . 



Inspection of the various terms in Eq. (|35p reveals that 
R is often much smaller than unity when the interaction 
is given by Eq. ([1]). Hence, a possibly convenient approx- 



imation to the coefficient Ac is from Eq. ([12 

B2 



Ac = Co + CiJ = Co + — y-^^Hs 



- 1 



(37) 



where J may be evaluated numerically or approximately 
with the help of Eqs. ([5))- piI)) . 

The full expression for the dynamical coefficient is writ- 
ten as 



Ar 



Co + ciJ + R, 



(38) 



using Eqs. ([TT|) and ([51]) . Now R from Eq. 
pressed as 



R w -0.055 J + 0.4745ii: - 9 x 10"'* {c^'^-^ - l) . (39) 

in view of Eqs. ([T|) and ([3]). Here we have introduced the 
function K for which we derive a convenient approxima- 
tion. 



K = 



12 
Ki 



Axx^ {x-2) (l-e-^^'- 
- (e^- -l)K2. 



where 



and 



Ki = 



Ax x^ 



{x~2)[ 



1-e" 



-Udh{x) 



Ko = 



2+<5 



dxx^ {x — 2)e 



(40) 



(41) 



(42) 



In the same spirit as in Ref. |32| . we approximate 
X (1 - e~^""(^)) w 2^0-^(^-2) - 2a^2e-2/.(a;-2)^ ^-^^-^ 

-. Wc then have 



C(A^ + l)(4-aO 
^3 



(43) 



where we have neglected the small term a£^^ /2fi'^. In the 
case of lysozyme at pH 4.5, the deviation of Eq. (|43|) from 
the exact result is smaller than about 3% for / > 0.05 
M and smaller than about 1% for / > 0.3 M. For 
the second integral we use the trapezoid approximation 



c2+5 

axg[X) RJ ^0 [gy. 
gleet a factor (1 + (5/2)^ 



J^^ dx g{x) w [5(2) + g{2 + 5)] {6 < 1) and we ne- 



K2 ~ 2(5^ cxp 



1 + 512 



(44) 



For lysozyme at pH 4.5 with 5 = 0.079 (see below), this 
approximation deviates less than about 5% from the ex- 
act value for / > 0.05 M and less than about 3% for 
/ > 0.2 M. 
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III. Comparison with experiment 

We compare our predictions of Ac as a function of the 
ionic strength I with experimental results for lysozymc at 
room temperature and at a pH of about 4.5. The added 
salt is NaCl and in most cases a small amount of Na ac- 
etate has been added as buffer. The reason for choosing 
lysozyme under these conditions is that we have previ- 
ously evaluated the range and strength of the short-range 
attraction^ and a lot of experimental data on the col- 
lective diffusion coefficient are available in the literature 
(see Fig. 131). 




□ 

-30 I ' ' ' 1 

0.5 1 1.5 2 

/(M) 

FIG. 3: Experimental data and theoretical predictions of Ac 
for lysozyme as a function of the ionic strength / at a pH of 
about 4.5. Data: black squares: Nystrom et al.— , pH 4.0, 
25 °C; grey squares: Mirarefi et al.— , pH 4.6; white squares: 
Mirarefi et al.— , pH 4.6; black diamonds: Muschol et al.— , 
pH 4.7, 20 °C; grey diamonds: Zhang et al.i^, pH 4.5, 20 
°C; white diamonds: Skouri et al.i^, pH 4.6, 20 °C; black 
triangles: Eberstein et al.— , pH 4.2, 20 °C; grey triangles: 
Leggio et al.— , pH 4.75, 25 °C; white triangles: Price et al.— , 
pH 4.6, 25 °C; black circles: Annunziata et al.— , pH 4.5, 25 
°C ; grey circles: Annunziata et al.— , pH 4.5, 25 °C; white 
circles: Retailleau et ali^, pH 4.0. In all cases, the supporting 
electrolyte is NaCl, often with a small amount of Na acetate 
added. The grey line denotes the theoretical curve setting 
i? = i.e. Eq. (|36|l with r given by Eq. dHJ , and the black 
line is the curve given by Eq. (|38|l . The functions J and K 
have been approximated as outlined in the text. 

Lysozyme has a moderate aspect ratio of about 1.5 
and we approximate it by a sphere of radius a = 1.7 
nm^. The dimensionless parameter ^ is then given by 
/i = 5.58%//, where the ionic strength / is given in M, 
and ^ = 0.209(Z/(1 -|-/.t))^. Here we follow our discussion 
in Ref. [H and use the adjusted charge on the lysozymc 
sphere Z = Z^/ / — 1 instead of the effective charge Zf.f / . 
Values of Z, Zf,ff and Z as & function of ionic strength 



can be found in Table I as well as the corresponding quan- 
tities fi and ^. For the range 5 and strength Ua of the 
attraction we use 6 = 0.079 and Ua = 3.70 which were 
computed on the basis of a wide variety of data on the 
second virial coefficient^. 

We next employ three methods to predict Ac theoret- 
ically. In the first, we compute r by equating the respec- 
tive second virial coefficients of Section II. B (sec Eqs. ([5]) 
and P^ ). Wc then calculate Ac from Eq. using 
Co — 1.454 and ci = 1.688. In the second method we use 
Eq. to determine Ac, where R is evaluated with the 
help of Eq. ([5^ . In both cases the approximations for 
J and K given by Eqs. ©-(HUl) and Eqs. dM]), (gSl) and 
(|44)) were used (see Table I and Fig. [S]). Note that there 
are no free parameters so the curves in Fig. [3] are pre- 
dictions not fits. For comparison, wc also calculate Ac 
from Eq. (jl4p exactly, that is by performing the integrals 
in Eqs. (fT5|) - p9)) numerically with the help of a highly 
accurate interpolation formula for h{x) (see Table I). Fi- 
nally, in Fig. Owe have also plotted data of Ac measured 
by several experimental groups. 



IV. Discussion 

In Section II. E we have outlined two approximate 
methods to calculate Ac- As one can see from Table 
I, both the direct method incorporating an approxima- 
tion for the residual R and the method relying solely on 
the stickiness r via the second virial yield results that are 
often close to the exact numerical computations. The di- 
rect method is, of course, somewhat more accurate. The 
T method breaks down below 0.2 M. Note that in the im- 
portant regime / > 0.2 M pertaining to protein crystal- 
lization, R is much smaller than the absolute magnitude 
of Ac- This may explain why Ac is a useful parameter 
to characterize the onset of crystallizationii. 

In Fig. [31 it is clear that there is a large degree of scat- 
ter which may be attributed to the systematic variation 
in sets of data from the various groups, especially at large 
ionic strengths (/ > 0.4 M). We do not know what is the 
cause of this. In one experimeni^i^, we do observe there is 
considerable scatter in a plot of the diffusion coefficient 
versus the protein solubility which might explain the ex- 
treme downturn of several data in Fig. [3] at about 0.5 
M. Fig. [3] also shows that our predicted curves lie fairly 
neatly in the midst of the swarm of data. We empha- 
size again that we have no adjustable parameters in our 
calculations except for a slight downward adjustment of 
the effective charge (sec also the discussion in Refs. [12 
and [ssh . The model is thus not inconsistent with the 
experimental data though we will have to await more ex- 
periments under conditions which are better controlled 
before one may reach a more definitive conclusion. In a 
similar vein, it is not possible to claim that the neglect 
of electrolyte friction assumed here is entirely warranted. 

In summary, wc have approximated proteins by spher- 
ical particles interacting by a hard-core and electrostatic 
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0.514 


0.286 


-2.48 


-1.97 


-1.96 


0.25 


10.2 


9.7 


8.7 


1.10 


2.79 


0.408 


0.193 


-4.38 


-3.97 


-3.95 


0.30 


10.2 


9.8 


8.8 


0.984 


3.06 


0.354 


0.155 


-5.81 


-5.46 


-5.45 


0.45 


10.3 


10.0 


9.0 


0.752 


3.74 


0.302 


0.109 


-8.82 


-8.53 


-8.51 


1.0 


10.4 


10.2 


9.2 


0.409 


5.58 


0.323 


0.0734 


-13.87 


-13.55 


-13.58 


1.5 


10.4 


10.3 


9.3 


0.295 


6.83 


0.349 


0.0655 


-15.72 


-15.38 


-15.43 


2.0 


10.4 


10.3 


9.3 


0.229 


7.89 


0.367 


0.0616 


-16.82 


-16.46 


-16.52 



TABLE I: Values of the actual charge Z of hen-egg-white lysozyme (from Ref . Isih . the effective charge Zefj (see Eq. Q), the 
lowered effective charge Z = Z^fi — 1, and dimensionless interaction parameters ^ and /i as a function of the ionic strength 
I. The pH equals 4.5 and ^ has been calculated using the lowered effective charge Z. R has been calculated from Eq. p9|l . 
T from Eq. (|12|l . Ac (via r) from Eq. p6|) and Ac (direct) from Eq. (|38}. In aU cases approximations for J and K given by 
Eqs. (|5)l- ([TU)) and (PO)). were used. The computation of the numerically exact Ac is explained in the text. 



repulsion together with a short-range attraction. An 
analysis of the two-particle statistics and hydrodynam- 
ics leads to a reasonable prediction of the ionic-strength 
dependence of the linear coefficient Ac- At high ionic 



strengths, when B2 is negative, the residual R is rela- 
tively small so there is then an interesting direct rela- 
tionship between Ac and B2 (Eq. ((57|) ) which could be 
tested experimentally. 
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